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SECTION  I 
SUMMARY 


The  problem  of  the  unsteady  transonic  flow  for  a  swept  wing  of 
moderate  sweep  is  formulated  in  the  small  disturbance  limit.  The  steady 
terms  in  the  resulting  potential  equation  contain  the  second  order  sweep 
terms  developed  at  NLR  (Amsterdam)  to  capture  the  important  swept 
shocks.  The  boundary  conditions  at  the  wing,  the  trailing  vortex  sheet, 
and  the  outer  computational  boundaries  are  given.  Ventilated  wall 
conditions  are  reviewed  as  well  as  the  resonant  wall  interference 
effects.  The  shortcomings  of  the  small  disturbance  hypothesis  is 
described,  suggesting  means  to  compensate  for  them.  Finally  the 
important  viscous  interactions  are  described  together  with  procedures  to 
incorporate  their  effects. 
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SECTION  II 
INTRODUCTION 


The  transonic  regime  Is  a  potential  troublespot  for  aeroelastlc 
problems  of  high  performance  aircraft.  This  Is  not  surprising  because  of 
the  well  known  sensitivity  of  the  flow  to  small  perturbations  In  the 
boundary  condition  at  the  wing  or  in  the  freestream.  The  consequence  of 
this  sensitivity  Is  magnified  by  the  presence  of  shock  waves,  the 
movement  of  which  can  cause  large  changes  of  the  lift  and  moment  on  the 
wing. 


Our  primary  objective  in  the  following  will  be  the  formulation  of 
the  unsteady  transonic  initial-boundary  value  problem  in  the  small 
disturbance  limit  applicable  for  aeroelastic  problems  of  swept  wings.  An 
essential  requisite  here  is  the  derivation  of  the  flow  equation  in  the 
proper  conservation  form  to  capture  the  important  moving  shock  waves. 
Such  an  objective  can  be  carried  out  in  a  straightforward  manner  using 
the  well  established  results  for  the  steady  case.  (See  References  1  and 
2.)  A  secondary,  but  still  important,  objective  is  to  review  the 
important  viscous  interactions  playing  an  essential  role  in  the 
positioning  and  movement  of  shock  waves,  and  to  describe  possible  means 
to  Incorporate  their  effects  in  a  manner  that  will  ultimately  not  add 
significantly  to  the  computer  time. 

In  a  formal  derivation  of  the  small  disturbance  problem,  a 
perturbation  potential  is  introduced  in  both  the  flow  equation  and  the 
boundary  conditions,  and  an  expansion  then  carried  out.  A  scaling  of  the 
independent  and  dependent  variables  and  a  physically  meaningful  limiting 
process  are  then  introduced  to  seek  a  nontrivial  problem  to  the  lowest 
order.  In  the  planar  steady  case,  the  result  was  the  familiar  problem 
embodying  the  classical  transonic  similarity.  The  same  similarity 
limiting  procedure  applied  to  the  three  dimensional  (3D)  case,  yielded  a 
3D  flow  equation  containing  only  the  iyyterm  to  represent  the  spanwise 
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effects  (Reference  3).  A  subsequent  examination  of  the  shock  jump 
condition  implied  by  this  equation  showed  that  only  shocks  essentially 
normal  to  the  frees tream  were  properly  captured,  precluding  the  proper 
capture  of  swept  shocks  (Reference  1).  By  adding  two  higher  order 
"sweep"  terms,  Bailey,  Ballhaus,  and  Lomax  (Reference  1)  showed  that  the 
proper  shock  condition  was  obtained  for  swept  shocks  arising  in  the  case 
of  a  constant  chord  yawed  wing  of  infinite  aspect  ratio.  Later  van  der 
Vooren,  Slooff,  Huizing,  and  van  Essen  (Reference  2)  by  a  different 
expansion  found  that  by  adding  the  same  sweep  terms  but  with  different 
(constant)  coefficients,  in  addition,  the  more  highly  swept  forward  shock 
typically  appearing  in  swept  wings  was  also  correctly  captured. 

The  above  higher  order  potential  equations  derived  in  References  1 
and  2  for  the  steady  3D  case  cannot  be  formally  derived  in  a  consistent 
way  by  a  similarity  limiting  process.  Thus  in  the  derivation  of  the 
unsteady  3D  equation  applicable  for  swept  wings,  we  shall  similarly 
dispense  with  limiting  processes  and  simply  add  the  necessary  second 
order  terms  to  capture  the  important  swept  shocks. 

In  the  following,  after  deriving  the  proper  flow  equation,  we  shall 
describe  the  boundary  condition  on  the  wing,  the  contact  conditions  on 
the  trailing  vortex  sheet,  and  the  outer  far-field  conditions.  Also  the 
ventilated  wall  condition  will  be  reviewed.  Limitations  of  the  resulting 
inviscid  small  disturbance  formulation  will  next  be  discussed,  with 
suggestions  for  their  compensation.  Finally,  the  dominant  viscous 
interactions  will  be  reviewed,  followed  by  means  to  incorporate  their 
effects  in  the  computations. 
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SECTION  III 

FORMULATION  OF  THE  PROBLEM 


In  the  present  section  we  shall  formulate  the  small  disturbance 
problem  starting  with  the  exact  potential  equation.  A  more  general 
formulation  is  given  in  more  detail  in  the  Appendix. 

1.0  BASIC  FLOW  EQUATIONS  (REFERENCES  a  &  5) 

Our  starting  point  is  the  exact  unsteady  potential  equation  given  by 

f?4+  $Wy+  =  0  (,} 


where 


Sea 


Here  subscripts  denote  partial  differentiations,  and 
(0  is  the  density  (slugs/ft3j 

the  velocity  potential  [ft2/secj 

X.Y2  Cartesian  coordinates  (X-streamwise;  Y-"spanwise" )  [ft] 
T  time  [sec] 

T  ratio  of  specific  heats 

tL  freestream  velocity  [it/s ec) 

M*,  freestream  Mach  number,  and 
freestream  density. 


Eq.  (1)  is  the  mass  conservation  equation  in  divergence  form,  and  it 
is  of  hyperbolic  type  independent  of  the  local  Mach  number. 

We  next  introduce  the  perturbation  potential  and  simultaneously 
nondimensionalize  the  variables  as  follows: 

$(Z,Y,Z. V  =  U»c  [x  *  *(x,y ,*,t] 

„  3  (3) 

21=  cx  Y=cy 

T-  (oc/ujt 

where  C  is  the  wing  reference  chord  [ft]  ,  and 

a  a  nondimensional  scaling  factor  on  the  time* 

The  perturbation  potential  is  next  introduced  into  Eqs.  (1)  and  (2), 
and  an  expansion  carried  out  retaining  only  those  quadratic  terms  needed 
to  yield  the  proper  steady  shock  jump  conditions  for  relevent  swept 
shocks.  We  obtain. 


V  •  -f  =  0 


where 


f  =  if,  +  Jf2  +  ^3 


-<5  (<5<?±+  %) 


f2  =  (•-*£%)  q>y 


(4) 


i3  =  %  H? 


* 

A  scaling  C  of  the  time  is  introduced  here  solely  to  identify  limiting 
cases.  For  the  computer  program  we  will  set  O  =  1. 
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h,T  ,3  ,k  are  the  unit  vectors  in  the  t,  x,  y  and  2  directions. 

Eq.  (4)  has  been  expressed  in  the  vector  form  to  facilitate  the 
derivation  of  the  shock  jump  conditions.  It  is  to  be  noted  that  this 
equation  ceases  to  be  valid  when  the  unsteadiness  becomes  too  severe 
(large  c)  as  in  the  case  of  a  “gust”  problem  modeled  by  an  abrupt  increase 
in  angle  of  attack.  Here  the  neglected  time  derivatives  as  the 
given  in  the  Appendix  will  have  a  significant  influence.  We  shall  see 
later  that  the  small  disturbance  wing  boundary  condition  will  also  become 
invalid  for  large <T. 

2.0  Shock  Jump  Conditions  (Reference  6) 

The  shock  surface  will  be  defined  by  s(x,y,z,t)  =  0.  The  normal 
vector  to  this  surface  will  then  be  given  by, 

_  ± 

n5=  (hst+  jsy  -*■  kse)  ( sf-tsf  +  sf)* 

On  s  =  0  we  have 

ds  =  dt  ■*  «SX  dx  +  Sy  dy  ■+  SjdZ  =  0  /$\ 


The  shock  velocity  components 


{UitVstWs)  are  given  by 


dz\ 

*tl„0 


so  that  the  above  equation  becomes 

st+  UsSx*  =  0  .  (6) 

The  shock  jump  conditions  are  based  upon  the  continuity  equation 

v-  f  =  0 
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together  with  the  supplementary  differentiability  conditions 


- 


£' 


- 


(%)t  =  (%K 

(%\  =  (%>y 

(%\  =  , 

or  in  "tensor"  notation 

V  M  =  0 

where 

M 

To  derive  the  shock  jump  conditions,  we  first  define  a  control 
surface  formed  by  the  upstream  and  downstream  faces  of  a  shock  surface 
element  and  introduce  Green's  formulae  equating  the  volume  integral 
of  7*  M  to  the  surface  integral  of  the  flux  of  M.  We  obtain 
directly 

LM]  •  ns=  0 


I 


—  _ 

where  [m]  *  M,  -  the  subscripts  1  and  2  denoting 

respectively  the  upstream  and  downstream  faces  of  the  shock;  that  is. 


[f,]sx  [fJ-Sy  ♦  ~  0 


Because  of  the  second  equation  in  Eq.  (7),  we  can  transfer  the  9-t 
occurring  in  "f  in  the  first  equation  to  the  term.  This  is  usually 
done,  as  in  LTRAN2  (Reference  31),  to  simplify  the  numerical  algorithm. 
Eqs.  (6)  and  (7)  then  constitute  the  shock  jump  conditions. 

It  will  be  reassuring  to  consider  next  several  familiar  simplified 
cases  -  namely,  the  one-dimensional  normal  shock  propagating  in  a 
uniform  steady  upstream  flow  and  the  steady  three-dimensional  (3D)  swept 
shock  case. 


In  the  case  of  the  unsteady  normal  shock,  Eqs.  (6)  and  (7)  simplify 


to 


S|  -+  Us  sx  =•  0 


Hi 


-5  ‘S'.Jst  +  =  0 

Kl  Sf  =  0 


Using  the  first  equation  in  the  second  and  third  equations  to 
eliminate  the  shape  derivatives  and  introducing  a  gal i lean 
transformation  £s.x-Us-t  with  we  obtain  the  jump 

condition 


*t,+  «V,(2-rU)} 


with 


Consider  next  the  steady  3D  case.  From  Eqs.  (5)  and  (7)  we  obtain 
the  jump  conditions 


8 


mm 


V  L%hs  0 


s 


(8) 


These  steady  jump  conditions  were  first  derived  by  van  der  Vooren, 
Slooff,  Huizing,  and  van  Essen  (Reference  2)  who  showed  that  the  jump 
condition  given  above  closely  approximated  the  exact  oblique  shock 
condition  for  both|(  the  yawed  wing  shock  as  well  as  for  the  more  highly 
swept  (and  thereby  weaker)  forward  shock.  (See  the  curves  labeled  as  NLR 
in  Figure  1  from  Reference  2). 

Bailey,  Ballhaus,  and  Lomax  (Reference  1)  earlier  derived  an 
alternate  steady  small  disturbance  equation  starting  with  the  exact  flow 
equations  in  a  different  form  than  Eq.  (1).  In  Reference  2  the  shock 
jump  conditions  implied  by  the  Bailey-Ballhaus  equation  were  also 
investigated  for  the  yawed  and  forward  shocks.  The  primary  difference 
with  respect  to  the  NLR  equation  was  a  less  accurate  capture  of  the 
forward  shock.  (See  curves  labeled  as  NASA  in  Figure  1.)  Note  also  in 
Figure  1  the  inadequate  jump  condition  for  the  classical  3D  equation  with 
only  the  term  (curve  labeled  Guderley-von  Karman). 

There  was,  however,  considerable  difficulty  in  capturing  the  weak 
forward  shock  in  a  steady  finite  difference  calculation  even  using  the 
more  exact  NLR  equation.  This  was  due  to  the  fact  that  such  highly  swept 
shocks  acquired  a  profile  spreading  over  many  more  mesh  points  than  the 
three  or  so  mesh  points  for  less  swept  stronger  shocks  as  a  result  of  the 
more  pronounced  diffusive  effects.  As  a  result,  even  using  a  relatively 
refined  mesh  about  the  forward  shock,  the  capture  was  still  not  fully 
satisfactory  as  can  be  seen  in  Figure  2. 
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On  the  other  hand,  there  is  usually  little  need  to  capture  the 
forward  shock  sharply,  since  the  important  primary  terminating  shock  is 
reasonably  well  captured  even  in  the  presence  of  a  smeared  and  poorly 
captured  forward  shock  *  This  can  be  seen  in  Figure  2  for  the  NIR 
equation  or  in  Figure  3  for  the  NASA  equation  for  the  same  ONERA  wing  at 
a  smaller  angle  of  attack. 

Thus  in  an  unsteady  calculation  the  improved  capture  of  the  forward 
shock  probably  is  not  warranted,  so  that  the  steady  terms  from  either  the 
NLR  or  the  NASA  equation  should  suffice.  But  since  the  NLR  equation  is 
no  more  complex  than  the  NASA  equation,  we  shall  adopt  the  steady  terms 
from  the  NLR  equation.  Thus  for  our  unsteady  equation,  we  shall  use 


=  a 

where  the  in  the  second  term  has  been  moved  to  the  first  term. 

It  is  clearly  evident  that  no  stretching  of  the  coordinates  can 
equalize  the  order  of  the  terms  in  this  equation.  More  significantly, 
similarity  laws  will  not  formally  exist  for  the  unsteady  3D  small 
disturbance  flows  described  by  Eq.  (9).  (See  Appendix.) 

In  the  case  of  the  purely  subcritical  or  purely  supersonic  flow,  the 
higher  order  terms  in  the  above  equation  are  indeed  higher  order  in  the 
absence  of  transonic  shocks,  and  one  obtains  in  the  limit  the  classical 
Prandtl-Glauert  or  the  Ackeret  linear  equation. 


*  Except  of  course  when  the  forward  and  rear  shocks  approach  each  other. 


Finally,  in  the  limit  consistent  with  the  small  disturbance 
approximation,  the  pressure  coefficient  is  given  by; 


T>-  P» 
iP.t'j 


-  2  (<px-+  <y  <Pt ) 


(10) 


3.0  COORDINATE  SYSTEM 

Consider  next  the  choice  of  the  coordinate  system  to  be  used  for  a 
general  swept  wing  in  a  finite  difference  procedure.  Here  the  experience 
gained  in  the  steady  case  will  be  directly  relevent.  Of  particular 
concern  here  is  to  select  a  coordinate  system  which  will  not  require  an 
excessive  number  of  mesh  planes  in  the  spanwise  direction.  Thus,  for 
example,  if  we  assume  a  cartesian  system,  it  is  clear  that  for  a  swept 
wing  the  spanwise  row  of  mesh  points  cutting  across  the  leading  edge  or 
across  a  shock  will  encounter  severe  spanwise  flow  gradients  as  in  the^ 
and  terms.  This  then  will  necessitate  a  fine  spanwise  spacing  of  the 
mesh  points.  Thus  in  the  steady  case,  Bailey  and  Ballhaus  (Reference  7) 
introduced  the  shearing  transformation  given  by 


*  ~ 
Cfy) 


*1  =  y 
4  =  * 


where  x  *  X|.t(y)  tf,e  equation  of  the  leading  edge  of  the  wing,  and 
c(y)  is  the  spanwise  chord  distribution.  The  above  transformation  maps 
the  swept  tapered  wing  in  the  (x,y,z)  space  to  a  rectangular  wing  in  the 
tt.*}. 0  space.  In  the  transformed  plane  the  mesh  configuration  is  now 
conformal  to  both  the  leading  and  trailing  edges,  so  that  a  reasonable 
mesh  spacing  In  the  spanwise  direction  can  be  used  with  the  resulting 
more  gradual  flow  variations  in  the  spanwise  coordinate. 
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Setting  <ff£  .oj  W*  ,y,2,-fj  •  Eq.  (9)  in  the  transformed 

plane  assumes  the  form  (see  Appendix  B) 

(cfj^  -*■  (C?3  =  0 


{  i 


-f.=  -<^Vc  V 


(11) 


Another  less  reliable  approach  with  respect  to  the  coordinate  system 
was  used  by  Schmidt,  Rohlfs,  and  Vanino  (Reference  8),  which  may  warrant 
a  reexamination  for  the  unsteady  case  due  to  its  simplicity.  Here  a 
cartesian  coordinate  system  is  used,  but  the  mesh  spacing  in  the  y 
direction  is  adjusted  such  that  the  relative  position  of  the  leading  edge 
relative  to  the  local  mesh  is  kept  uniform  over  the  span  to  moderate  the 
consequences  of  using  an  inadequate  mesh  in  the  spanwise  direction.  Here 
the  exact  solution  can  be  used  in  the  steady  case  to  tailor  the  mesh 
pattern.  The  advantage  here  will  be  the  possibility  of  using  the  simpler 
Eq.  (9)  rather  than  the  above  Eq.  (1^1) . 


4.0  WING  BOUNDARY  CONDITION 

It  is  well  known  that  in  boundary  value  problems,  the  boundary 
conditions  invariably  play  a  key  role.  This  is  particularly  the  case  for 
the  unsteady  transonic  problem  where  the  flow  tangency  condition  on  the 
wing  serves  as  the  driving  or  forcing  function  for  the  flow  development. 
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At  the  surface  of  the  wing  defined  by  B(R,T)  =  0  where  R  is  the 
space  vector  R(X,Y,Z),  the  normal  component  of  local  flow  velocity  vector  ,  t 

( )  relative  to  the  moving  surface  element  with  velocity  Vg  must  be 
zero.  That  Is, 


(V4>-Ve)-  N,  =  o. 


(12) 


Here  is  the  instantaneous  outward  normal  to  the  wing  surface,  B(R,T)-0, 
given  by 


where  the  gradient  V  operates  only  with  R 


Further,  by  differentiation,  we  have 


Replacing  Vj  by  7<p  in  this  equation  using  Eq 
required  boundary  condition  for  the  wing;  that  is 


we  obtain  the 


If  we  now  introduce  the  perturbation  potential  and  the 
nondimensional  variables  from  Eq.  (3)  into  Eq.  (13)  and  retain  only  the 
lowest  order  terms  in  the  perturbation  potential,  we  obtain 

%(x.y  >',t)  =  o’  bt°(x,y,t )  •+ 
where  B(X,Y,Z,T)  =  Z- 

c* 

Expanding  (p2  in  a  Taylor  series  about  z  =  0,  we  obtain 


<f>2i(*-y.o,t)b%  •••  =  a  K  (i4) 


or  to  the  lowest  order 


(15) 


Eq.  (15)  constitutes  the  important  wing  boundary  condition. 

It  should  be  noted  here  that  for  wings  with  blunt  leading  edges,  not 
only  is  the  assumption  of  small  perturbations  invalid,  but  the  neglected 
higher  order  terms  of  the  Taylor  expansion  on  the  left  side  of  Eq.  (14) 
invariably  overwhelm  the  retained  lowest  order  term.  Moreover,  omission 
of  these  higher  order  terms  also  renders  pure  plunging  motion  degenerate 
for  ->>  0.  Finally,  note  that  the  small  disturbance  approximation  is 
violated  if  the  time  scaling  O'  is  taken  too  large. 

Consider  next  the  important  problem  of  the  gust.  Classically  such  a 
case  was  treated  for  example  by  simply  imposing  a  moving  step  change  of 


the  angle  of  attack.  No  doubt  a  measure  of  the  gust  loading  can  be 
obtained  in  this  way,  but  it  would  be  of  importance  to  treat  the  gust  in 
a  more  realistic  manner,  at  the  least,  to  assess  the  classical  gust 
treatment.  Thus  for  this  purpose,  a  gust  packet  with  a  “typical"  profile 
of  upwash  in  both  space  and  time  is  introduced  into  the  flow  at  the  far 
upstream  boundary  and  allowed  to  convect  downstream  toward  the  wing. 
Here  some  experimentation  is  required  to  compensate  for  the  diffusion  of 
the  gust  packet,  such  that  the  desired  gust  is  achieved  just  upstream  of 
the  wing.  Such  experimentation  should  be  logically  carried  out  in  the 
planar  case.  Such  an  approach  would  parallel  that  carried  out 
experimentally  by  M.  Platzer  of  the  Navy  Postgraduate  School  who  employed 
an  unsteady  jet  flap  upstream  of  the  test  model  to  generate  a  gust. 

Finally,  one  must  impose  the  important  Kutta  condition  at  the  wing 
trailing  edge  which  controls  the  local  chordwise  circulation  and 
therefore  the  local  shed  rate  of  vorticity.  That  is,  we  have  the 
requirement  that  as  the  trailing  edge  is  approached  along  the  upper  and 
lower  surfaces  of  the  wing,  the  pressure  approaches  a  common  value.  Thus 
we  have  the  condition 

[cf]  =  Cp(*re,y ,0+,i)~  cf(xn,y.o:t)*  o  (w) 

where  x  =  x^  (y)  defines  the  trailing  edge  and  z  =  0+,  O'represent 
respectively  the  upper  and  lower  surface  of  the  wing;  or  using  Eq. (10), 
we  obtain  for  the  Kutta  condition, 

(h+ch)M=  o.  (i7, 


5.0  TRAILING  WAKE  CONDITION 

The  physical  picture  as  given  by  Helmholtz's  vorticity  law  is  that 
at  a  given  point  on  the  trailing  edge,  as  the  local  chordwise  circulation 
changes,  the  change  in  circulation  must  be  shed  into  the  flow  at  the 
trailing  edge.  Such  shed  vorticity  is  transferred  to  the  fluid  particle 
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in  residence  at  the  trailing  edge  at  the  time  of  the  shedding,  thereafter 
traveling  downstream  always  tied  to  the  fluid  particle. 

In  an  Eulerian  representation,  along  the  trailing  vortex  sheet,  one 
must  have  a  continuity  of  both  streamline  slope  and  pressure  across  the 
sheet.  These  conditions  will  be  imposed,  as  in  the  wing  tangency 
condition,  in  a  quasi-planar  fashion  on  the  z  =  0  plane.  That  is, 

Oz]  =  9afx,y,0tt)-  <P2(x,y,0~t)  =  0  (18) 

and 

[C„]  *  Cp(xtyJo+t)-Ct>(*,y)0->i)=  0  . 

The  second  condition  above  can  be  reexpressed  in  terms  of  the  perturbation 
potential  using  Eq.  (10)  as  follows: 

(h) 

or  in  terms  of  the  general  integral 

[<p(*j y,-fcj  =  . 

The  substantive  or  convective  derivative  in  Eq.  (19)  indicates  that 
the  jump  in  generated  at  the  trailing  edge  travels  downstream  along  the 
particle  path  line  at  the  particle  speed  1  /o  as  required  by  the  Helmholtz 
vortex  law.  (Note  here  that  the  apparent  distortion  of  the  particle  speed 
is  due  to  the  change  of  scale  of  the  time  introduced  in  Eq.  (3).)  Thus 
the  W  to  be  assigned  at  a  point  x  -  xTS  =  x0  downstream  of  the  trailing 
edge  at  t»i#  Is  that  which  was  shed  at  the  trailing  edge  at  the  earlier 
time  t  *  t„  -  x.ff. 


As  G+>0  (1  o v/  unsteadiness  or  in  the  low  frequency  limit  in  the  case 
of  harmonically  oscillating  flow),  we  have  l<p(x  ,  y,t)J  =  L9>(x7e  yjj. 

That  is,  l<pj  for  a  given  y  and  t  along  the  wake  is  invariant  with  x  at 

M  »l 

the  value  at  the  trailing  edge.  In  this  low  frequency  limit,  the  flow 
changes  so  slowly  that  in  the  time  required  for  the  shed  vorticity  to 
travel  far  downstream,  the  flow,  in  particular  the  t<pJ  at  the  wing  trailing 
edge,  has  not  significantly  changed. 

6.0  OUTER  BOUNDARY  CONDITIONS 

The  far  field  conditions  for  the  unsteady  3D  problem  are  that 
sufficiently  far  upstream  and  laterally  from  the  wing,  the  flow  must 
approach  the  freestream  conditions,  while  far  downstream  it  must  approach 
a  Trefftz  type  flow,  characterized  by  the  trailing  vortex  system.  In  a 
calculational  procedure,  outer  boundaries  at  finite  but  large  distances 
from  the  wing  must  be  used,  and  the  question  then  arises  as  to  what 

conditions  to  prescribe  on  such  computational  boundaries. 

To  answer  this  question,  let  us  first  recall  that  in  the  case  of 

unsteady  30  flow,  we  have  a  hyperbolic  problem.  Thus,  for  example,  the 
"Mach  cone"  for  the  low  frequency  planar  case  is  given  in  Figure  4  from 
Reference  9.  Guidance  as  to  the  necessary  far  field  condition  for  the 

unsteady  3D  case  can  thus  be  obtained  by  examining  the  simplified  case  of 
a  steady  supersonic  flow  over  an  axial  symmetric  body  shown  in  Figure  5. 
Here  the  lines  abc,  cb"d,  bef,  and  eb"  are  characteristics  or  Mach  waves. 
Thus  if  we  take  an  outer  boundary  no  closer  to  the  body  than  c'cc" ,  one 
can  determine  the  flow  at  the  body  independent  of  the  condition  prescribed 
along  cc".  We  can  for  example  set  the  flow  sloped,  =  0  along  the  boundary 
c'cc".  If,  however,  we  take  the  closer  boundary  b'bb",  then  the  flow 
along  fd  on  the  body  will  be  distorted  from  the  correct  flow  unless  the 
right  condition  is  prescribed  along  bb"  of  the  outer  boundary.  In  the 

past  a  condition  along  bb"  has  been  sought  assuming  some  form  of  an 
"outgoing  wave"  or  radiating  condition.  It'  is  clear  that  such  an  approach 
cannot  be  strictly  exact  since  the  flow  along  bb"  is  clearly  dependent  not 
only  on  the  initial  values  along  be,  but  on  the  characteristic  initial 
values  along  be,  beyond  the  outer  boundary  bb".  On  the  other  hand,  if  the 
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boundary  b'bb"  is  located  sufficiently  remote  from  the  body,  a  condition 
such  as  the  above  outgoing  wave  condition  or  the  simpler  condition  <P„=  0 
(solid  wall  condition)  can  be  prescribed  without  serious  errors  being 
reflected  back  onto  the  body  along  the  segment  fd. 

Thus  in  the  present  unsteady  3D  transonic  problem  we  shall 
prescribe  the  condition  <p„  =  0  on  the  lateral  boundaries  taken  sufficiently 
far  from  the  wing  that  either  the  reflected  perturbations  do  not  have 
sufficient  time  to  impinge  back  onto  the  wing,  or  if  they  do,  the 

reflected  errors  are  negligible.  This  can  be  easily  verified  in  the 

calculations  by  successively  extending  the  boundaries  outward  until  the 

desired  result  is  achieved.  Frequently  mesh  stretching  is  used  to  extend 
the  location  of  the  outer  boundaries.  Here  an  additional  source  of  error 
can  arise  if  the  stretched  mesh  is  no  longer  able  to  resolve  properly  the 
outgoing  unsteady  waves,  but  in  general  such  an  error  probably  will  not 
have  a  significant  effect  in  the  above  circumstances. 

The  far  downstream  boundary  condition  requires  special  consideration 
since  the  flow  in  the  far  downstream  boundary  is  affected  by  the  unsteady 
trailing  vortex  sheet  passing  through  it.  (See  Figure  6.)  If  the 

downstream  boundary  is  taken  sufficiently  far  from  the  wing,  the  flow  in 
the  downstream  boundary  will  be  independent  of  the  influence  of  the  wing 
and  will  be  a  function  only  of  the  trailing  vortex  sheet  fluxing  through 
it  with  frozen  vorticity  traveling  at  the  propagation  speed  l/c.  (See 
Section  5.0)  That  is,  the  flow  in  a  plane  transverse  to  the  freestream 
moving  at  the  velocity  1/c  will  be  invariant  with  time,  that  is  <p(x,y,z,t) 
=  •Ht-ffx.y),  or 


(20) 


This  then  is  the  condition  to  be  prescribed  on  the  far  downstream 
boundary,  and  it  is  equivalent  to  Cp  =  0.  In  the  limit  this 

condition  reduces  to  f*  =  0. 
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Finally  note  that  if  the  downstream  boundary  is  not  taken 
sufficiently  far  downstream,  the  effect  of  the  wing  cannot  be  ignored,  and 
the  flow  in  a  plane  transverse  to  the  freestream,  traveling  downstream  at 
a  velocity  of  1/(T,  will  not  be  invariant  with  time.  In  such  a  case,  the 
use  of  a  boundary  condition  such  as  Eq.  (20)  in  'a’cComputation  could  lead 
to  numerical  errors  or  instabilities  if  the  proper  outfluxing  of  the 
l  trailing  vorticity  is  not  allowed. 

7.0  VENTILATED  WIND  TUNNEL  WALL  CONDITIONS 

In  the  present  section  we  shall  review  the  wall  boundary  conditions 
for  a  ventilated  wind  tunnel  with  slotted  or  perforated  walls.  The 
conditions  for  the  slotted  wall  are  given,  for  example,  in  Reference  10; 
while  those  for  the  perforated  walls  are  given  in  Reference  11.  These 
conditions  understandably  are  highly  idealized,  and  their  validity  has  yet 
to  be  demonstrated  in  the  steady  case  in  a  real  wind  tunnel  environment. 
These  conditions  will  be  extrapolated  to  the  unsteady  case  by  assuming 
their  applicability  to  the  instantaneous  flows. 

7.1  Slotted  Wall 

The  idealized  condition  for  the  slotted  wall  is  based  on  the 
assumption  that  the  pressure  at  the  slot  is  equal  to  the  plenum  pressure, 
while  the  normal  velocity  is  zero  along  the  solid  portion  of  the  wall. 
These  mixed  conditions  are  converted  to  a  homogenized  condition  applicable 
on  a  substitute  boundary  just  interior  to  the  actual  walls.  This 
condition  is  derived  by  applying  the  conservation  of  mass  between  the 
actual  and  substitute  boundaries  assuming  the  streamwise  variation  of  the 
flow  to  be  negligible.  In  the  usual  case  that  the  spacing  d  between  the 
slots  is  large  relative  to  the  slot  width,  and  with  the  plenum  pressure 
set  equal  to  the  freestream  pressure,  one  obtains  the  slotted  wall 
condition 


Here  CP)W  and  8W  are  the  pressure  coefficient  and  the  streamline  slope  on 
the  substitute  boundary,  which  can  be  taken  at  the  actual  wall  location; 
and  K  is  a  constant  dependent  only  on  the  slot  configuration.  In  the  3D 
case  0w  is  the  slope  of  the  streamline  in  a  streamwise  plane  normal  to  the 
wall.  The  primary  source  of  error  in  the  derivation  of  the  above 
condition  is  the  assumption  of  inviscid  flow. 

The  above  condition  is  valid  only  along  those  portions  of  the  wall 
where  an  outflow  from  the  tunnel  to  the  plenum  occurs.  When  an  inflow 
occurs,  the  total  pressure  of  the  incoming  air  has  suffered  a  drop  due  to 
the  mixing  of  the  moving  tunnel  air  with  the  stagnated  plenum  air.  In 
this  case  in  Reference  10  an  alternate  condition  was  derived,  again  with 
simplifying  assumptions,  given  by 

Cp)W "  1C  -  Cf), plenum  (22) 


Here  Cp  p/e„um  is  the  plenum  pressure  coefficient  assumed  as  constant,  and 
K’  is  the  essentially  constant  perturbation  potential  at  the  slot. 

In  the  unsteady  case  we  assume  that  Eqs.  (21)  and  (22)  apply  to  the 
instantaneous  flow,  setting 

C*  =  *  2  a  <?t)  . 


The  outstanding  problem  with  these  ideal  slotted  wall  conditions  is 
the  selection  of  the  constants  K  and  K'  since  in  a  lifting  case,  both 
outflow  as  well  as  inflow  can  be  expected.  Moreover,  one  faces  the 
physically  unreasonable  discontinuity  of  the  boundary  condition  as  the 
wall  flow  switches  between  an  inflow  and  an  outflow. 

Understandably,  in  steady  wall  interference  studies,  the  outflow 
condition  Eq.  (21)  Is  usually  assumed  to  prevail  along  the  entire  wall,  so 
that  only  the  tunnel  dependent  K  is  needed. 


Theoretical  as  well  as  experimental  derivations  of  K  have  been 
reviewed  by  Barnwell  (Reference  10)  who  also  recommended  values  of  K  as  a 
function  of  the  ratio  of  the  width  to  spacing  of  the  slots. 

7.2  Perforated  Wall 

The  basis  of  the  perforated  wall  condition  is  that  the  mass  flux 
through  the  wall  is  proportional  to  the  pressure  difference  across  the 
wall.  This  then  leads  to  the  perforated  wall  condition 

Cf>tW=- 2(%ta<f>t)  =  ,  (23) 


Here  the  plenum  pressure  has  been  set  equal  to  the  freestream  pressure, 
and  the  constant  R  is  a  measure  of  the  wall  porosity.  The  parameter  R  is 
usually  determined  in  a  laboratory  setup  by  measuring  the  mass  flux 
through  the  perforated  plate  as  a  function  of  the  pressure  drop  across  the 
plate  without  an  ambient  flow  along  the  pl«te.  How  applicable  such  a 
determination  of  R  is  for  the  case  of  the  wind  tunnel  wall  is  not  clear, 
since  in  the  latter,  one  has  a  complex  boundary  layer  flow  passing 
transversely  across  the  perforations. 

As  for  the  slotted  wall,  Eq.  (23)  no  longer  applies  in  the  case  of 
inflow  from  the  plenum  to  the  tunnel.  Here  again  the  inflow  has  a  greatly 
reduced  total  pressure  (approximately  equal  to  the  plenum  pressure),  and 
this  low  energy  air  then  forms  a  buffer  layer  along  the  wall  that  softens 
the  wall  effect.  Some  consequences  of  this  are  discussed  in  Reference  12. 

Clearly  the  effects  of  this  low  energy  buffer  layer  are  difficult  to 
assess,  so  that  the  outflow  condition  (23)  is  usually  assumed  to  apply 
over  the  entire  wall.  Mokry  (Reference  11)  found  that  improved  prediction 
of  the  wall  interference  is  obtained  if  a  different  value  of  R  is  used  for 
the  upper  and  lower  walls,  somewhat  compensating  for  the  less  accurate 
treatment  of  the  inflow. 


7.3  Summary  of  the  Wall  Boundary  Conditions 

There  is  no  question  that  the  ventilated  wall  conditions  given  above 
are  idealized  conditions  that  hopefully  yield  ttie  "first  order"  effects  of 
the  wall  interference.  In  summary  we  have  therefore  the  following  wall 
conditions. 

Slotted  Wall: 


+  K%3  =  0 

where  K  is  the  slot  parameter  given,  for  example,  in  Figure  4  of  Reference 

10. 


Perforated  Wall: 

where  typical  values  of  the  porosity  R  are  given  in  Reference  13. 

For  completeness  we  have 
Solid  Wall: 


Free  Jet: 

%+9%  ~  0 

Side  walls  of  the  wind  tunnel  are  usually  closed  walls,  so  that  the 
condition  iy*o  Is  prescribed  on  these  walls. 

7.4  Wind  Tunnel  Wall  Interference 

Assuming  that  the  ventilated  wall  conditions  given  above  are  valid. 
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wall  interference  is  obtained  by  comparing  the  computational  results  with 
the  wall  condition  and  with  an  infinite  freestream  condition.  If 
experimental  pressure  distribution  data  are  available  in  the  wind  tunnel 
of  interest,  they  can  be  used  to  determine  the  wall  parameters. 

In  using  experimental  results  in  this  fashion,  one  must  check  that 
the  test  section  length,  both  upstream  and  downstream  of  the  wing,  is 
sufficiently  large  that  the  effects  of  the  wing  have  subsided  sufficiently 
that  freestream  conditions  prevail  at  the  extremities  of  the  test 
section.  In  particular  an  erroneous  determination  of  the  freestream  Mach 
number  can  arise  if  the  upstream  effects  of  the  wing  have  not  sufficiently 
attenuated  where  the  measurement  for  the  freestream  velocity  is  taken. 

Finally,  we  shall  review  the  resonant  wall  interference  that  arises 
in  unsteady  testing  which  can  significantly  distort,  for  example,  the 
results  in  flutter  testing.  Here  the  unsteady  waves  generated  by  the  wing 
resonate  between  the  wing  and  the  wind  tunnel  wall  at  discrete  frequencies 
that  are  functions  of  the  velocity  of  sound  and  the  ratio  of  the  tunnel 
height  to  wing  chord.  A  typical  wind  tunnel  result  showing  these 
interference  effects  from  Reference  14  is  shown  in  Figure  7  for  the  case 
of  the  pitching  oscillation  of  the  NACA  65-010  airfoil.  Here  it  is  seen 
at  the  resonant  frequency  that  the  amplitude  of  the  lift  and  moment 
variations  drops  abruptly,  while  a  large  phase  lag  increase  is  seen 
carrying  over  an  extended  range  of  frequencies  below  resonance.  The 
results  of  Reference  14  further  showed  the  resonant  frequency  fell  in  the 
low  frequency  range  that  is  relevent  to  the  flutter  phenomena. 

The  above  results  are  for  a  planar  wing  in  subcritical  flow  with 
solid  wind  tunnel  walls.  Thus  in  the  unsteady  testing  of  a  swept  wing  in 
the  supercritical  range  with  ventilated  wind  tunnel  walls,  the  possibility 
of  such  resonant  wall  interference  must  be  kept  in  mind.  Of  particular 
importance  is  the  possible  large  distortion  of  the  phase  which  may  be 
present  well  below  the  resonant  frequency  where  the  distortion  of  the  lift 
and  moment  amplitudes  is  not  evident.  Here  the  calculational  procedure 
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based  on  the  unsteady  transonic  3D  problem  formulated  in  the  previous 
sections  would  be  invaluable  to  check  and  correct  for  the  above  resonant 
wall  interference.  Viscous  interaction  effects,  to  be  described  in 
Section  10,  must  be  incorporated  into  these  calculations. 

8.0  SUMMARY  OF  THE  PROBLEM  FORMULATION 

The  purpose  of  the  present  section  is  to  summarize  the 
initial-boundary  value  problem  in  the  small  disturbance  limit  derived  in 
detail  in  the  previous  sections  and  to  define  several  limiting  cases. 

Thus  the  general  flow  equation  for  the  perturbation  potential  is 


-<J  (oy  2  %  )t  +  (fi\-  A%2-  ig.; 


=  0 


where 


I- 
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The  significance  of  the  underlining  will  be  discussed  shortly.  The 
pressure  coefficient  is 

Cp=  -2(%  +  <T9t). 

•  •  •  •  • 

The  boundary  condition  on  the  wing  (z  =  t>(x,y,t))  is 

<P2(x,yO*t)  =  b,Vx,y,t;+  <5  b^(x,yt±) 


Additionally  the  Kutta  condition 
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is  imposed  at  the  subsonic  trailing  edges. 


Along  the  trailing  vortex  sheet  we  have  the  conditions 

i 

02]=  o 

0 

Finally  for  the  outer  conditions  on  boundaries  taken  sufficiently  far 
from  the  wing,  we  have. 


^=0  (far  upstream  and  far  laterally) 

^dx+C|t)<P=  0  (far  downstream) 

To  complete  the  formulation,  an  initial  starting  flow  must  be 
prescribed,  which  can  be  taken  as  any  readily  available  steady  flow  or  the 
uniform  flow  <J>  =  0  if  such  a  flow  is  unavailable. 

There  are  several  limiting  cases  of  the  above  problem.  In  the  limit 
that  O’  *-0,  we  have  the  fluid  mechanically  defined  quasi-steady  flow  in 
which  terms  underlined  by  the  solid  line  in  the  above  equations  are 
deleted.  Here  the  time  appears  only  parametrically  in  the  wing  condition. 

There  Is  also  the  quasi-steady  problem  as  defined  by  the 
aeroelastician.  Here  the  time  derivatives  in  the  flow  equation  are 
deleted,  but  they  are  retained  in  all  the  other  equations. 

We  can  next  define  the  low  unsteadiness  problem  in  which  the  terms 
underlined  by  the  dashed  lines  are  deleted.  Here  the  unsteadiness  is 
manifested  only  through  the  term  in  the  potential  equation,  and 

implicitly  in  a  parametric  fashion  through  the  wing  boundary  condition. 


Finally,  we  can  define  a  quasi-high  unsteadiness  problem  in  which 
only  the  term  is  deleted  in  the  potential  equation  with  all  time 
derivatives  retained  in  the  other  equations. 

What  are  now  the  implications  and  significances  of  the  above  various 
limiting  cases?  For  the  general  case  including  the  gust  problem,  no  terms 
in  the  above  formulation  can  be  deleted.  Moreover,  solutions  to  the 
complete  problem  are  essential  to  establish  the  range  of  validity  of  the 
approximate  cases  defined  above,  in  particular  to  establish  upper  bounds 
on  the  permissible  values  of  6" . 

The  quasi-high  unsteadiness  case  is  of  particular  interest.  Here  the 
addition  of  the  time  derivatives  in  the  boundary  conditions  can  be  easily 
implemented  in  the  low  unsteadiness  case,  to  then  extend  the  range  of 
applicability  to  cover  most  of  the  unsteady  cases  of  interest  except  for 
the  gust  problem.  (Preliminary  results  of  van  der  Vooren  of  the  NLR. ) 

It  is  appropriate  finally  to  comment  on  the  simplified  approach  not 
considered  above  in  which  the  unsteadiness  is  considered  as  a  small 
perturbation  about  a  steady  background  flow.  Such  a  linearized  unsteady 
approach  is  undoubtedly  meaningful  so  long  as  significant  shock  waves  or 
more  precisely  their  significant  motion  do  not  arise.  Clearly,  if  one 
considers  the  region  of  the  flow  including  the  wing  surface  traversed  by 
the  shocks,  the  unsteadiness  here  cannot  be  considered  a  small 
perturbation  since  the  flow  changes  abruptly  from  a  supersonic  to  subsonic 
flow  (or  vice  versa)  as  the  shock  travels  over  points  of  such  a  subdomain. 

It  may  be  recalled  that  the  expected  transonic  aeroelastic  problems 
are  aggravated  by  the  magnified  lift  and  moment  changes  due  to  the  motion 
of  shock  waves,  most  probably  to  the  phase  lags  between  the  shock  motion 
and  the  wing  motion.  There  simply  Is  no  possibility  to  capture  the  shock 
motion  correctly  with  the  above  time  linearized  approach. 
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In  a  similar  vein,  there  have  been  serious  queries  regarding  the 
possibility  that  the  flow  for  a  wing  in  combined  plunging  and  pitching 
oscillatory  motions  can  be  considered  as  a  superposition  of  the  flow  with 
a  pure  plunging  motion  and  a  pure  pitching  oscillation.  Again  if 
significant  shock  motion  arises,  the  flow  is  essentially  nonlinear, 
precluding  such  a  superposition. 


SECTION  IV 

LIMITATIONS  OF  AND  COMPENSATION  FOR  THE  SMALL  DISTURBANCE  APPROXIMATION 


The  nature  of  the  errors  caused  by  the  small  disturbance 
approximation  in  the  unsteady  case  can  be  directly  extrapolated  from  the 
steady  case  which  has  been  extensively  studied.  Thus  at  blunt  leading 
edges,  not  only  is  the  small  disturbance  potential  equation  inaccurate, 
but  the  tangency  condition  is  grossly  in  error  as  well  as  the 
quasi-planar  fashion  in  which  it  is  imposed.  Here  the  important 
consequence  is,  not  the  inaccuracy  about  the  stagnation  point,  but  the 
accuracy  of  the  subsequent  upper  surface  leading  edge  expansion  to 
supersonic  velocities.  This  establishes  the  suction  plateau  and  hence 
the  strength  of  the  shock  waves.  Elsewhere  on  the  wing  where  the  surface 
slopes  become  excessive  as  at  large  values  of  the  angle  of  attack  or  in 
the  aft  camber  region  of  supercritical  airfoils  at  low  or  moderate  angle 
of  attack,  significant  errors  can  also  arise. 

In  the  unsteady  case  of  a  swept  wing,  there  is  a  simple  means  to 
compensate  for  these  inaccuracies  (Reference  15).  Here  the  steady  exact 
flow  over  a  suitably  selected  representative  planar  airfoil  or  the  3D 
wing  itself  can  be  used  to  fashion  shape  changes  in  the  nose  or  aft 
camber  regions  of  the  wing  to  compensate  for  these  shortcomings. 

Thus  consider  a  swept  wing  with  a  moderate  sweep  (~30°)  with  a 
supercritical  airfoil  section.  First  a  characteristic  planar  flow  is 
defined  at  the  midsemispan  using  local  sweep  theory  and  assuming  uniform 
downwash.  This  2D  flow  is  then  computed  by  an  existing  exact  potential 
method.  The  2D  small  disturbance  code  Is  then  used,  distorting  the  nose 
slopes  in  a  smooth  fashion  as  with  the  Rlegels  rule  (Reference  16)  and 
adding  aft  decambering  ramp,  such  that  the  exact  pressures  are  obtained. 
If  the  2D  small  disturbance  code  has  a  capability  to  prescribe  pressures 
In  lieu  of  the  surface  slopes,  the  latter  aft  decambering  ramps  can  be 
evolved  directly  by  prescribing  the  exact  pressures  In  the  aft  portion  of 
the  upper  and  lower  surfaces. 


The  airfoil  shape  changes  thus  evolved  are  then  suitably 
incorporated  into  the  swept  wing.  Clearly  some  development  is  required 
in  the  technique,  using  the  exact  3D  flow  as  a  validation  base.  In  the 
case  of  a  rigid  body  pitching  oscillation  case,  the  above  shape  tailoring 
would  he  carried  out  at  the  mean  angle  of  attack. 

It  is  clearly  not  the  intent  of  the  above  procedure  to  obtain  a 
"perfect"  compensation  for  the  instantaneous  pressure  distributions.  It 
is  desirable,  however,  to  obtain  an  acceptable  prediction  of  the  shock 
motion  and  the  resulting  lift  and  moment  time  histories  acceptable  for 
aeroelastic  applications. 
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SECTION  V 

TURBULENT  VISCOUS  INTERACTIONS 


In  the  present  section  we  shall  review  the  important  viscous 
interactions  for  the  unsteady  transonic  flow  over  wings  of  moderate 
sweep.  Since  there  is  little  data  available  for  the  unsteady  case,  we 
must  therefore  extrapolate  our  experiences  with  steady  flows.  We  shall 
then  describe  two  methods  to  incorporate  these  viscous  interactions  into 
the  inviscid  flow  methods.  They  are  the  viscous  ramp  method  (Reference 
17)  and  the  Lag  Entrainment  method  (References  18  and  19),  which  have 

been  selected  for  their  ability  to  predict  the  important  viscous  effects 
without  intolerably  complicating  the  overall  procedure. 

For  swept  wings  of  moderate  sweep  (£30o)  and  aspect  ratios  in  the 
range  of  five  or  greater,  we  shall  postulate  that  the  viscous 
interactions  can  be  treated  in  a  suitable  quasi-planar  fashion  in  cuts 
normal  to  the  sweep.  That  this  is  a  reasonable  approximation  can  be  seen 
in  Figure  8  for  the  case  of  a  yawed  wing  of  constant  chord  and  infinite 
aspect  ratio  of  35°  sweep  tested  by  Elsenaar,  van  den  Berg,  and 

Lindhout  (Reference  20).  Here  it  is  seen  that  the  displacement  thickness 
computed  in  the  above  quasi-planar  way  closely  approximates  the  3D 

calculations  with  McLean's  method  (Reference  21).  Also  seen  here  is  the 
inappropriateness  of  taking  cuts  in  the  streamwise  direction.  Similar 
results  were  obtained  by  Nash  and  Tseng  (Reference  22).  The  above 
approximation  should  be  acceptable  over  most  of  the  wing,  except  in  the 
tip  region,  when  severe  separation  is  absent.  With  severe  separation, 
particularly  near  the  trailing  edge,  a  spanwise  flow  of  the  "dead  air" 
occurs  greatly  increasing  the  viscous  displacement  in  the  tipward 
direction,  thereby  invalidating  the  quasi-planar  model. 

In  planar  unsteady  flows,  it  has  been  made  plausible  in  Reference  17 
that  the  important  viscous  interactions  for  sufficiently  moderate 
unsteadiness  can  be  treated  in  a  quasi-steady  way  in  a  suitable 

framework.  (We  shall  review  this  in  more  detail  shortly.)  Thus  a 
discription  of  the  viscous  interactions  in  the  steady  planar  case  will  be 
relevant. 


In  Figure  9  we  compare  schematically  the  flows  over  a  supercritical 
airfoil  and  the  corresponding  pressure  distributions  for  the  case  of 
inviscid  and  viscous  flows.  The  dominant  and  well  known  effect  of 
viscosity  shown  here  is  the  reduction  of  the  shock  pressure  rise  and  the 
consequent  upstream  displacement  of  the  shock  in  the  viscous  case 
relative  to  the  inviscid  case.  Here  the  reduction  of  the  viscous  shock 
pressure  rise  is  due  to  the  abrupt  increase  of  the  boundary  layer 
displacement  (a  wedging  effect)  just  downstream  of  the  shock  which 
converts  the  normal  shock  in  the  inviscid  case  to  the  observed  weaker 
oblique  shock  with  its  reduced  pressure  rise.  The  wedging  displacement 
at  the  shock  is  the  consequence  of  the  severe  reduction  of  the  fullness 
of  the  boundary  layer  velocity  profile  and  increase  of  the  boundary  layer 
thickness  after  passage  through  the  abrupt  shock  adverse  pressure 
gradient.  The  reduced  shock  pressure  rise  then  displaces  the  shock 
upstream. 

A  second  effect  of  viscosity  also  shown  here  is  the  reduction  of  the 
plateau  loading  caused  by  the  aft  decambering  effect  of  the  boundary 
layer  and  near  wake.  This  effect  is  particularly  significant  in  the  case 
of  aft  cambered  airfoils  since  the  resulting  aft  decambering  has  a  global 
influence  on  the  pressure  distribution  as  shown  above. 

In  summary,  the  primary  viscous  effect  is  the  shock-boundary  layer 
interaction.  For  aft  cambered  airfoils,  additionally  the  aft 
displacement  effects  of  the  boundary  layer  and  near  wake  will  be 
important.  An  acceptable  determination  of  the  lift  and  moment  on  the 
wing  cannot  be  made  without  incorporating  these  effects. 

Two  methods  to  calculate  the  above  viscous  effects  for  the  unsteady 
3D  case  are  the  viscous  ramp  method  and  the  lag  entrainment  method. 
These  methods  are  essentially  planar  steady  methods.  They  will  be 
relevent  for  the  3D  unsteady  case  since,  as  described  earlier,  a  quasi  2D 
flow  Is  postulated  and  the  resulting  planar  flows  then  considered 
quasi -steady  in  the  proper  framework.  These  methods  have  been  selected 
for  their  simplicity  for  near-term  applicability. 
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1.0  The  Viscous  Ramp  Method 

This  procedure  is  a  phenomenological  method  in  which  a  priori 
determined  shape  changes  simulating  the  viscous  displacement  effects  are 
incorporated  into  the  inviscid  procedure.  Here  the  shock-boundary  layer 
interaction  is  modeled  by  placing  a  wedge-nosed  ramp  at  the  base  of  the 
shock  to  obtain  the  reduced  shock  pressure  rise.  The  appropriate  wedge 
angle  is  determined  via  the  shock  polar  using  a  postulated  post-shock 
pressure  based  upon  measurements.  (See  Figure  10.)  Other  features  of 
the  wedge  ramp  as  well  as  its  placement  relative  to  the  shock  profile  are 
developed  by  calculations  using  as  a  comparison  base  suitable  steady 
measured  pressure  distributions.  Models  of  such  ramps  for  example  have 
been  developed  by  Murman  (Reference  23)  who  used  the  shock  detachment 
value  of  the  wedge  angle,  and  in  Reference  17. 

In  the  unsteady  case,  the  shock  in  general  will  be  out  of  phase  with 
the  motion  of  the  airfoil.  Thus  the  wedge  ramp  must  be  tied,  not  to  the 
instantaneous  angle  of  attack,  but  to  the  prevailing  shock,  with  the 
wedge  angle  adjusting  instantaneously  to  the  instantaneous  shock 
strength.  In  other  words  the  boundary  layer  just  downstream  of  the 
shock,  and  hence  the  wedge  ramp,  is  assumed  to  equilibrate  in  a  time 
sufficiently  short  that  the  shock  does  not  change  significantly  in 
position  or  strength.  This  approximation  is  not  unreasonable  provided 
the  flow  unsteadiness  is  not  excessive.  Thus  consider  the  pitching 
oscillation  case  with  the  reduced  frequency  k  =  0.468  treated  in 
Reference  17.  Here  the  characteristic  time  for  the  unsteadiness  of  the 
flow,  and  hence  that  of  the  shock  motion,  is  the  period  of  the 
oscillation,  which  in  the  above  case  is  13.4  t:.,  where  x,  is  the  time 
required  for  the  airfoil  to  travel  one  chord  at  the  freestream  velocity. 
On  the  other  hand  the  characteristic  time  for  the  equilibration  of  the 
boundary  layer  just  downstream  of  the  shock  to  the  changing  shock  should 
be  no  greater  than  onex„.  Thus  we  can  tailor  the  wedge  ramp  to  the 
instantaneous  shock  in  a  quasi-steady  fashion. 

Consider  next  the  aft  displacement  effects  of  the  boundary  layer  and 
the  near  wake.  These  effects  are  less  dramatic  than  those  of  the 


shock-boundary  layer  interaction,  but  nevertheless  they  still  can  play  an 
important  role,  even  affecting  the  shock  motion  by  modifying  the 
pressures  upstream  of  the  shock.  In  an  unsteady  case,  the  aft  pressure 
gradients  determining  the  boundary  layer  displacement  are  closely  in 
phase  with  the  instantaneous  angle  of  attack.  We  shall  therefore 
incorporate  the  aft  viscous  ramps  in  a  quasi-steady  fashion  to  the 
instantaneous  angle  of  attack. 

The  shaping  of  the  aft  ramps  will  be  carried  out  in  a  steady  flow  at 
a  characteristic  angle  of  attack  for  the  unsteady  flow.  Thus  in  the  case 
of  the  pitching  oscillation,  the  mean  angle  of  attack  will  be  used.  For 
this  purpose  the  steady  experimental  pressure  distribution  for  the 
airfoil  at  the  appropriate  Mach  number  and  Reynolds  number  is  required. 
The  shaping  of  the  aft  ramps  is  then  determined  by  an  inviscid 
calculation  in  which  the  measured  pressures  are  prescribed  as  boundary 
condition  in  lieu  of  the  surface  slopes.  The  results  will  yield  aft 
shape  changes  required  to  yield  the  measured  pressures,  the  latter 
incorporating  the  effects  of  the  boundary  layer  as  well  as  the  near 
wake.  The  resulting  aft  ramps  are  then  frozen  onto  the  airfoil.  They 
will  then  be  off-design  at  angles  of  attack  other  than  the  mean  value, 
but  the  errors  will  be  assumed  tolerable  since  variations  of  the  aft 
pressure  gradients,  or  more  precisely  their  effect  on  the  boundary  layer 
displacement,  are  usually  not  large. 

The  above  decoupled  treatment  of  the  viscous  interactions  at  the 
shock  and  in  the  aft  region  is  permitted  only  when  the  Reynolds  number  is 
sufficiently  high  that  the  type  B  interaction  defined  by  Pearcey 
(Reference  24)  does  not  arise.  That  is,  the  boundary  layer  is 
sufficiently  thin  that  it  recovers  from  its  encounter  with  the  shock 
adverse  pressure  gradient  prior  to  its  encounter  with  the  second  aft 
adverse  subcritical  pressure  recovery.  Without  a  full  recovery  of  the 
loss  of  the  velocity  profile  fullness  through  the  shock  before 
encountering  the  second  adverse  gradient,  the  effects  of  the  latter  are 
greatly  magnified,  precluding  the  assumed  decoupled  treatment  of  the 
separate  interactions. 


In  addition,  if  the  boundary  layer  thickness  upstream  of  the  shock 
is  "unnaturally"  large,  for  example,  as  the  result  of  an  overly 

conservative  tripping  of  the  boundary  layer  in  a  wind  tunnel  test,  a  type 
B  interaction  can  also  arise  if  the  increase  of  the  thickness  of  the 

recovered  boundary  layer  is  sufficiently  large. 

The  above  viscous  ramp  procedure  consequently  reduces  to  the 

computation  of  an  inviscid  flow  with  a  priori  determinable  time  dependent 
shape  changes.  In  the  case  of  the  wedge-nosed  ramp  for  the 

shock-boundary  layer  interaction,  when  the  Mach  number  upstream  of  the 
shock  is  1.25,  the  wedge  angle  corresponds,  in  the  case  of  the  shock 
detachment  angle,  to  approximately  5.3o.  The  movement  of  such  a  wedge 
ramp  on  the  surface  of  the  airfoil  creates  a  significant  disturbance  to 
the  highly  sensitive  flow.  Thus  in  the  implementation  of  the  viscous 
ramp  procedure  into  a  numerical  procedure,  skill  must  be  exercised  to 
avoid  numerical  instabilities  and  excessive  additional  computing  time. 

2.0  Lag  Entrainment  Method 

The  viscous  ramp  method  described  above  is  a  phenomenological 
procedure,  devoid  of  differential  equations,  requiring  minimal  additional 
computing  time  for  near-term  use.  A  more  refined  procedure  will  require 
appropriate  flow  differential  equations  to  describe  the  viscous  flow,  and 
such  equations  must  be  able  to  treat  shock-boundary  layer  interactions 
and  locally  separated  flows  in  a  reasonable  fashion.  The  turbulent 
transport  model  applicable  for  such  flows  are  presently  not  available, 
and  prospects  for  their  development  are  discouraging.  We  must  therefore 
select  the  level  of  description  of  the  viscous  flow  such  that  this 
shortcoming  can  be  most  simply  circumvented  in  a  phenomenological 
manner.  We  must  further  continually  keep  in  mind  the  requirement  of 
minimal  computing  times.  Since  we  shall  consider  the  3D  unsteady  flow  to 
be  quasi  2D,  and  the  latter  quasi-steady,  we  need  consider  only  the 
steady  2D  case.  These  considerations  then  lead  naturally  to  integral 
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methods  with  ordinary  differential  equations,  rather  than  partial 
differential  equation  methods;  in  particular  to  the  lag  entrainment 
methods  which  have  had  success  in  previous  steady  2D  cases  (see 
References  18,  19  and  25). 

In  treating  the  viscous  interactions  with  the  integral  equation 
method,  one  considers  the  8  ,  the  viscous  displacement,  a  quantity 
required  to  input  directly  the  viscous  effects  into  the  inviscid 
procedure.  However,  for  the  concept  of  the  displacement  thickness  to  be 
strictly  meaningful,  the  pressure  must  be  invariant  across  the  viscous 
layer.  Such  an  assumption  is  most  probably  violated  in  the  region  of  the 
shock  impingement  and  about  the  trailing  edge  where  the  boundary  layer 
passes  into  the  wake.  To  circumvent  this  conceptual  difficulty,  we  shall 
define  a  fictitious  boundary  layer  for  which  the  concept  of  S*  is 
meaningful,  tailoring  the  necessary  empirical  inputs  for  the  turbulent 
transport  to  fit  the  framework  of  this  substitute  boundary  layer. 

The  starting  point  for  the  lag  entrainment  method  is  the  integral 
form  of  the  streamwise  momentum  and  the  continuity  equations  obtained  by 
an  integration  across  the  boundary  layer;  that  is, 

Streamwise  momentum: 
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Continuity: 
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(25) 


Here  x  is  the  streamwise  coordinate;  Ue  the  streamwise  velocity,  V.  the 
transverse  velocity,  and  the  density  all  at  the  outer  edge  of  the 
boundary  layer;  the  displacement  or  mass  flux  defect  thickness  is 
defined  by 

S 
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where  £  is  the  boundary  layer  thickness;  the  momentum  thickness  is 
defined  by 

* .  -i 

the  form  factor  H  =  S/6  ;  and  fy*  the  wall  skin  friction 

coefficient.  For  convenience,  the  above  lengths  are  nondimensionalized 
by  the  wing  chord,  the  velocities  by  the  freestream  velocity,  and  the 
density  by  its  freestream  value.  Note  also  that  the  velocities  are  not 
perturbations. 


Eqs.  (24)  and  (25)  contain  the  six  unknowns  £  ,  £  ,  9  ,  ,  Ve  and  Cf  . 

(Here  Me  and  pe  are  obtained  from  Ue  from  the  isentropic  relations.)  For 
a  fully  determined  set  of  equations,  four  additional  supplementary 
equations  are  required.  These  are  furnished  by  a  skin  friction  law,  by 
defining  a  two  parameter  velocity  profile  law,  and  by  a  lag  entrainment 
hypothesis  furnished  by  Bradshaw  (Reference  26).  The  final  required 
equation  is  furnished  from  the  interaction  with  the  outer  inviscid  flow. 
These  auxiliary  relations  are  next  discussed. 


For  the  skin  friction  law,  several  forms  have  proven  to  be 
satisfactory.  Thus,  for  example,  in  References  19  and  22,  the 
Spaulding-Chi  law  was  used,  while  in  Reference  18,  the  Ludwieg-Til Imann 
formulae  has  been  used.  These  laws  are  given  in  the  functional  form 

> 

Cf  =  Cf(Me,T?0,H)  (26) 

where  Rfl  is  the  local  Reynolds  number  based  upon  0,  the  momentum 
thickness. 

For  the  velocity  profiles  we  shall  assume  the  two  parameter  family 
developed  in  Reference  27,  which  was  shown  to  be  adequate  to  represent 
the  diverse  range  of  profiles  measured  by  Seddon  (Reference  28)  for  the 
case  of  the  normal  shock  impingement  on  the  turbulent  boundary  layer 
along  a  flat  plate.  Here  Cole's  law  of  the  wake  is  used  together  with 
van  Driest's  law  of  the  wall  with  a  laminar  sublayer.  Using  the 
Crocco-Busemann  law  for  zero  heat  transfer  to  express  the  density  profile 
in  terms  of  the  above  velocity  profile,  one  can  obtain  the  displacement 
and  momentum  thicknesses  in  terms  of  the  two  parameters  g  and  C^;  that  is, 

s=  skcf)  0=  8<i'cf)  <27> 

Ce  =  Ve  /Ue  in  Eq.  (25)  is  defined  as  the  entrainment  coefficient 
which  is  equal  to  the  slope  of  the  streamline  of  the  inviscid  flow  just 
entering  the  boundary  layer.  It  is  a  measure  of  the  rate  at  which  outer 
inviscid  flow  is  being  entrained  locally  into  the  boundary  layer.  Its 
introduction  into  the  boundary  layer  analysis  was  first  due  to  Head 
(Reference  29). 


For  the  skin  friction  law,  several  forms  have  proven  to  be 
satisfactory.  Thus,  for  example,  in  References  19  and  22,  the 
Spaul ding-Chi  law  was  used,  while  In  Reference  18,  the  Ludwieg-Til lmann 
formulae  has  been  used.  These  laws  are  given  in  the  functional  form 

Cf  =  Cf(>VVH)  (26) 


where  Rfl  Is  the  local  Reynolds  number  based  upon  0,  the  momentum 
thickness. 

For  the  velocity  profiles  we  shall  assume  the  two  parameter  family 
developed  in  Reference  27,  which  was  shown  to  be  adequate  to  represent 
the  diverse  range  of  profiles  measured  by  Seddon  (Reference  28)  for  the 
case  of  the  normal  shock  impingement  on  the  turbulent  boundary  layer 
along  a  flat  plate.  Here  Cole's  law  of  the  wake  is  used  together  with 
van  Driest's  law  of  the  wall  with  a  laminar  sublayer.  Using  the 
Crocco-Busemann  law  for  zero  heat  transfer  to  express  the  density  profile 
in  terms  of  the  above  velocity  profile,  one  can  obtain  the  displacement 
and  momentum  thicknesses  in  terms  of  the  two  parameters  $  and  C^;  that  is, 


$*=  S(S,Cf)  6=  B(S,Cf)  (27) 


C£  =  Ve  /U4  in  Eq.  (25)  is  defined  as  the  entrainment  coefficient 
which  is  equal  to  the  slope  of  the  streamline  of  the  inviscid  flow  just 
entering  the  boundary  layer.  It  is  a  measure  of  the  rate  at  which  outer 
Inviscid  flow  is  being  entrained  locally  into  the  boundary  layer.  Its 
Introduction  into  the  boundary  layer  analysis  was  first  due  to  Head 
(Reference  29). 
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Bradshaw  (Reference  26)  empirically  related  the  entrainment 
coefficient  functionally  with  the  maximum  shear  stress  within  the  local 
boundary  layer  profile  by  the  relation 


Ce 


10 


'HIM 


(28) 


Prior  history  effects  in  the  turbulent  entrainment  are  incorporated 
by  assuming  a  simple  lag  entrainment  equation,  which  in  Reference  18  is 
given  by 


—~x )  = 

A.  / 

T««\ 
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i  '  ft 

fe  > 

(29) 


Here  A.  is  a  lag  constant  to  be  empirically  determined,  and  ' l is 
the  maximum  shear  stress  for  the  equilibrium  boundary  layer  in  which  the 
velocity  and  shear  stress  profiles  are  invariant  with  x. 

To  determine  the  maximum  shear  stress  from  the  postulated  velocity 
profile,  in  Reference  18  an  eddy  diffusivity  is  assumed.  Though  the 
maximum  shear  stress  can  be  determined  once  the  eddy  diffusivity  and  the 
velocity  profile  are  assumed,  in  Reference  18  an  empirically  determined 
location  of  the  maximum  shear  stress  within  the  profile  is  used,  which  is 
different  for  attached  and  separated  flows.  The  resulting  expressions 
for  the  equilibrium  case  are  given  in  Figure  7  of  Reference  18. 

e* 

Finally  Ue  or  d  is  obtained  from  the  inviscid  flow.  When  Is 
furnished  by  the  inviscid  flow,  a  direct  problem  for  the  boundary  layer 
is  defined  for  which  the  separation  point,  Cf  =  0,  is  a  singular  point. 
When  the  £  is  inputed,  an  indirect  problem  is  defined,  but  now  the 
separation  point  is  no  longer  a  singular  point.  Quite  clearly  the  latter 
formulation  for  the  boundary  layer  is  preferred. 
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Thus  in  summary,  Eqs.  ( 24 ) - (29) ,  and  an  input  of  or  d£/dx  or  a 
from  the  inviscid  flow  calculations  constitute  the  system  of  equations 
for  the  boundary  layer,  in  particular  the  formulation  from  Reference  18. 
Green  in  Reference  19  proposed  a  similar  lag  entrainment  procedure  which 
differed  from  the  above  primarily  in  the  expressions  for  the  lag 
entrainment  and  in  the  choice  of  the  velocity  profile.  Empiricism  has 
been  incorporated  in  the  two  methods  in  different  ways,  so  that  a  direct 
comparison  cannot  be  made  without  comparing  the  results  for  specific 
cases. 

The  validity  of  Green's  lag  entrainment  method  has  been  examined  in 
Reference  25  for  a  number  of  cases  by  comparisons  with  experiments.  In 
particular  in  Figure  11  we  reproduce  one  of  the  shock-boundary  layer 
Interaction  cases  of  interest,  namely  that  tested  by  Seddon  (Reference 
28)  for  the  normal  shock  impingement  on  a  flat  plate.  Seddon  measured 
the  velocity  profiles  and  the  surface  pressures  and  from  these  determined 
all  of  the  dependent  terms  arising  in  Green's  equations.  In  the  results 
of  Figure  11,  Green's  method  was  used  to  calculate  Me,0,  and  Cj  inputing 
the  measured  $*.  The  computed  quantities  are  then  compared  with  Seddon 's 
experimental  results  in  Figure  10.  The  agreement  of  0  and  Cf  are  near 
perfect,  with  the  main  discrepancy  arising  in  the  post-shock  recovery  of 
Me.  Most  plausible  cause  for  the  latter  discrepancy  is  the  displacement 
effects  of  the  side  wall  boundary  layers  which  were  not  removed.  Similar 
good  agreements  with  experiments  were  obtained  for  other  cases. 

In  the  above  indirect  problem  for  the  boundary  layer,  where  S  is 
inputed  from  the  inviscid  flow,  it  has  been  demonstrated  that  Green's 
equations  can  be  Integrated  in  a  stable  fashion  through  the  shock 
pressure  rise  without  difficulty.  These  are  encouraging  results  that 
recommend  the  use  of  the  lag  entrainment  methods. 

3.0  Viscid-Inviscld  Flow  Coupling 

In  the  previous  sections  we  have  formulated  the  problems  for  the 
inviscid  flow  and  the  boundary  layer  flow.  We  must  next  consider  the 
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coupling  of  these  two  flows.  In  the  classical  low  speed  case,  this  was 
accomplished  by  the  well  known  weak  interaction  iterative  procedure.  The 
starting  point  for  the  iteration  is  an  inviscid  calculation  from  which 
the  surface  pressure  distribution  required  for  the  boundary  layer 
calculation  is  obtained.  The  latter  calculation  then  determines  the 
displacement  thickness  which,  when  added  to  the  original  airfoil  shape, 
results  in  a  new  effective  airfoil  shape  and  a  new  starting  point.  The 
above  procedure  is  repeated  until  an  invariant  pressure  distribution  is 
obtained. 

In  the  low  speed  case  the  convergence  of  the  iteration  is  usually 
rapid  so  long  as  significant  separation  does  not  arise.  With  separation, 
a  strong  interaction  arises  requiring  the  boundary  layer  and  the  inviscid 
flow  to  be  treated  simultaneously. 

The  viscous  interaction  in  the  case  of  the  supercritical  flow  over 
airfoils  is  an  extreme  case  of  a  strong  interaction.  However,  due  to  the 
extremeness  of  the  shock-boundary  layer  interaction  and  the  sensitivity 
of  the  transonic  inviscid  flow,  in  general,  to  perturbations  of  the 
surface  boundary  condition,  a  direct  coupling  alone  may  not  preclude 
stability  difficulties  in  the  resulting  numerical  procedure. 

The  problem  of  the  transonic  coupling  has  been  considered,  for 
example,  by  Le  Balleur  (Reference  30)  who  treated  the  outer  inviscid  flow 
with  a  steady  small  disturbance  relaxation  procedure  and  the  boundary 
layer  with  an  integral  method.  A  semi -inverse  iterative  procedure  with 
under-relaxation  was  employed  for  the  coupling  where  the  inviscid  flow 
was  treated  as  a  direct  problem  (surface  slope  prescribed  to  yield  the 
pressure)  and  the  boundary  layer  treated  as  an  indirect  problem  (  $  or 
slope  prescribed  to  yield  the  pressure). 

The  iterative  procedure  was  as  follows.  With  an  initial  slope 
distribution  given,  the  pressure  distribution  was  computed  for  both 
flows.  The  difference  of  the  resulting  pressure  distributions  was  then 
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used  as  a  forcing  function  to  evolve  a  slope  correction  by  an  approximate 
locally  linearized  procedure,  the  correction  then  suitably  under-relaxed 
to  serve  as  a  new  point  of  departure. 

The  movement  of  the  resulting  viscous  ramp  across  mesh  points  from 
iteration  to  iteration  in  the  above  direct  formulation  of  the  inviscid 
problem  undoubtedly  strains  the  stability  of  the  inviscid  flow 
algorithm.  This  was  relieved  by  Le  Balleur  who  embedded  a  highly  refined 
mesh  about  the  shock  which  moved  with  the  shock.  Thus  in  this  way  with  a 
moving  mesh  a  given  surface  mesh  point  is  spared  the  sudden  confrontation 
with  a  large  change  of  the  boundary  condition,  as  would  arise  in  the  case 
of  a  point  in  a  fixed  mesh  when  the  viscous  wedge  moved  onto  it. 

For  the  unsteady  problem,  the  procedure  of  Le  Balleur  has  the 
undesirable  feature  of  increased  computing  times  due  to  the  need  of  a 
highly  refined  moving  mesh  about  the  shock  and  the  need  to  iterate  at 
each  time  step  to  achieve  a  time-accurate  solution. 

One  has  the  further  deficiency  that  the  inviscid  flow  relaxation 

procedure  is  incapable  of  treating  the  shock  profile  region  due  to  the 

large  truncation  errors.  The  consequent  large  numerical  dissipative 

effects  are  primarily  responsible  for  the  shock  profile.  In  Le  Balleur' s 

£ 

result  near  the  surface  the  concavity  of  the  $  surface  widened  the  shock 
profile  over  many  more  mesh  points  than  the  customary  three  or  so  that 
would  usually  prevail  and  which  must  prevail  outside  the  region  of 
influence  of  the  concavity.  In  the  above  widened  portion  of  the  shock 
profile,  the  use  of  a  refined  mesh  has  undoubtedly  reduced  the  truncation 
errors,  but  one  must  then  question  the  validity  of  the  small  disturbance 
equation  in  this  region. 

In  the  following  we  shall  suggest  an  alternate  procedure  in  an 
attempt  to  circumvent  the  above  objections,  coupling  directly  the 
inviscid  and  boundary  layer  flows.  The  shock  profile  region  is  then 
treated  by  an  ad  hoc  procedure  in  which  a  bridging  condition  for  the 
boundary  layer  is  obtained  to  apply  across  the  shock  profile  region. 
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Consider  the  low  frequency  case  where  the  LTRAN  2  code  developed  by 
Ballhaus  and  Goorjian  (Reference  31)  may  be  used  for  the  inviscid  flow. 
This  method  employs  a  two  step  implicit  approximate  factorization 
difference  algorithm  in  which  the  first  step  is  a  predictor  step  with  the 
time  and  x  (streamwise)  derivatives  treated  implicitly  and  the  y 
derivative  treated  explicitly.  In  this  step  (the  x-sweep)  the  value  of 
the  potential  along  a  given  row  of  mesh  points  is  determined 
simultaneously  imposing  the  far  upstream  and  downstream  boundary 
conditions. 

In  the  second  step,  the  corrector  step,  corrections  to  the  potential 
obtained  in  the  x-sweep  are  obtained  by  treating  the  time  and  y 
derivatives  implicitly  in  y-sweeps,  imposing  timely  conditions  along  the 
airfoil . 

For  the  boundary  layer,  the  steady  lag  entrainment  method  will  be 
used  assuming  the  boundary  layer  to  be  quasi-steady  relative  to  the 
instantaneous  inviscid  flow. 

Formally  a  simultaneous  coupling  of  the  inviscid  and  boundary  layer 
flows  is  first  accomplished  in  the  following  way  using  Figure  12.  Here 
we  show  the  pressure  distribution  at  the  prior  time  (top),  the  inviscid 
mesh  system  (middle),  and  the  boundary  layer  flow  for  the  shock  region 
(bottom).  The  shock  profile  region  has  been  shaded. 

Thus,  in  the  unsteady  inviscid  flow,  assume  that  the  x-sweeps  at  the 
new  time  have  been  accomplished  as  well  as  the  subsequent  y-sweeps 
through  the  column  labeled  0.  The  boundary  layer  has  been  also 
determined  to  station  0.  We  now  consider  the  determination  of  the  flow 
at  station  1,  coupling  directly  the  inviscid  and  boundary  layer  flows. 

From  the  boundary  layer  difference  equations,  an  algebraic  relation 
between  the  velocity  Ue  (or  the  pressure)  and  d£/dx  at  the  station  1  is 
obtained  in  terms  of  the  known  conditions  at  station  0.  This  relation  is 
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now  to  apply  as  the  boundary  condition  at  the  mesh  point  1  for  the 

inviscid  flow.  Since  the  pressure  is  proportional  to  <f>K  in  the  low 
frequency  limit,  the  prescribed  %  along  the  airfoil  (y  =  0)  can  be 

integrated  to  yield  <p  at  the  station  1  in  terms  of  <j ?  at  the  stations  0 

and  1  to  second  order  accuracy.  Thus  the  boundary  layer  equations  yield 

a  boundary  condition  of  the  type  <p  +  afy  =  b  for  the  inviscid  flow;  that 
is  a  condition  of  the  third  type. 

In  an  actual  calculation  a  smaller  Ax  would  be  required  than  shown 
in  Figure  12  for  the  boundary  layer  calculations  within  the  shock  profile 
as  well  as  also  for  the  inviscid  flow  to  effect  the  direct  coupling. 
Thus  one  undesirability  of  Le  Balleur's  procedure  for  the  unsteady 
problem,  namely  the  need  of  a  fine  mesh  about  the  shock  profile,  is  not 
yet  avoided.  To  remedy  this,  as  well  as  to  seek  a  more  rational  and  less 
vulnerable  procedure  through  the  shock  profile,  an  ad  hoc  procedure  is 
suggested  for  the  shock  profile  region. 

It  may  be  recalled  that  in  the  case  of  inviscid  flow,  the  shock 
profile  region  was  excluded  from  the  physical  domain.  Thus  a  fine  mesh 
was  used  about  the  shock  to  minimize  this  loss  of  the  flow  domain.  Also, 
a  conservative  difference  algorithm ,  insuring  a  telescoping  of  the 
truncation  errors,  was  used  to  insure  the  correct  jump  condition  across 
corresponding  points  of  the  excluded  region. 

The  boundary  layer  underlying  the  shock  profile  will  be 
correspondingly  treated  in  a  consistent  manner,  the  goal  to  obtain 
bridging  conditions  across  the  questionable  region  under  the  shock 
profile.  Such  a  bridging  condition  would  relate  the  boundary  layer 
conditions  just  downstream  of  the  shock  to  the  known  state  just  upstream 
of  the  shock. 

Thus  consider  again  Figure  12.  Based  upon  the  pressure  distribution 
at  the  prior  time  (top  of  Figure  12),  the  pre-shock  and  post-shock  points 
are  defined  relative  to  the  shock  profile.  At  the  downstream  station  the 
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pressure  is  assumed  to  be  the  empirically  observed  value  as  a  function  of 
M,  ,  the  Mach  number  at  the  upstream  station,  as  used  in  the  viscous  ramp 
procedure.  The  shock  polar  is  then  utilized  to  obtain  the  corresponding 
wedge  angle,  for  example  the  shock  detachment  value.  The  d<f/dx  is  then 
assumed  to  vary  in  a  prescribed  manner  between  this  value  and  the  known 
value  at  the  upstream  station.  (Here  some  experimentation  with  the 
boundary  layer  code  will  be  required  to  evolve  a  distribution  of  S 
through  the  shock  profile  to  achieve  the  desired  post-shock  pressure.) 
With  the  resulting  &  prescribed,  the  well-behaved  indirect  boundary 
layer  problem  (see  Reference  25)  is  solved  to  yield  the  pressure 
distribution  through  the  shock  profile  as  well  as  the  boundary  layer 
conditions  at  the  post-shock  station.  The  determination  of  the  latter 
constitutes  the  bridging  condition  mentioned  earlier.  The  resulting 
velocity  distribution  converted  to  a  potential  distribution  is  then  used 
as  a  boundary  condition  for  the  inviscid  flow  y-sweeps  through  the  shock 
profile. 

The  continuation  of  the  y-sweeps  downstream  of  the  shock  is  treated 
in  the  coupled  manner  described  earlier  obtaining  the  starting  condition 
for  the  boundary  layer  from  the  bridging  condition. 

Here  to  forestall  possible  instabilities  far  downstream,  it  may  be 
necessary  to  replace  the  viscous  equations  by  patching  on  a  suitable 

asymptotic  similarity  wake  at  a  sufficiently  large  distance  (perhaps  two 
chords)  downstream  of  the  trailing  edge.  Further,  the  displacement 
thickness  of  such  an  asymptotic  wake  may  not  decrease  rapidly  enough  to 
be  compatible  with  the  freestream  conditions  prescribed  at  the  downstream 
computational  boundary.  A  distortion  of  the  asymptotic  wake  somewhat 
upstream  of  the  downstream  boundary  would  then  be  required  such  that 
d$*/dx  =  0  over  perhaps  a  distance  of  one  chord  upstream  of  the 

boundary.  These  alterations  of  the  wake  occur  sufficiently  far 

downstream  that  they  should  have  little  effect  on  the  surface  pressure 

distributions. 
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In  the  above  procedure,  in  bridging  the  shock  profile,  the  d  £/dx 
obtained  from  the  inviscid  flow  will  not  in  general  agree  with  the 
assumed  distribution  used  for  the  boundary  layer.  This  mismatch  might  be 
eliminated  by  an  iterative  procedure,  but  this  is  probably  not  warranted 
due  to  the  questionable  validity  of  the  boundary  layer  equations  which 
assumed  a  constancy  of  the  pressure  across  the  boundary  layer  as  well  as 
by  the  use  of  a  vulnerable  turbulent  entrainment.  The  inviscid  flow  is 
also  inadequately  treated  in  this  region  as  mentioned  earlier. 

The  procedure  as  described  above  will  not  require  a  refined  mesh 
about  the  shock  in  the  inviscid  flow  nor  iterations  at  each  time  step. 
The  procedure  described  above  is  a  reasonable  approach,  but  its  viability 
must  await  numerical  verification  and  possible  further  development.  Its 
formulation  has  served  as  a  vehicle  to  focus  upon  the  nature  of  the 
difficulties  inherent  in  the  problem. 
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SECTION  VI 
SUMMARY  REMARKS 


It  is  essential  that  the  unsteady  fluid  dynamic  code  for  use  in 
important  aeroelastic  problems  of  swept  wings  be  capable  of  capturing  the 
important  unsteady  shock  waves.  To  accomplish  this,  the  appropriate 
equation  in  proper  conservation  form  must  be  used,  and  the  significant 
viscous  interactions  at  the  shock  and  in  the  aft  region  of  the  wing  must 
be  incorporated.  It  is  desirable  to  incorporate  the  viscous  effects  with 
minimum  additional  complexity. 

Thus  time-linearized  approximation  to  the  flow  equation  will  be 
inadequate  for  important  supercritical  flows. 

Following  the  steady  state  developments,  second  order  sweep  terms 
have  been  added  to  the  classical  3D  small  disturbance  potential  equation 
to  insure  the  correct  shock  jump  conditions.  Such  an  equation  does  not 
admit  similarity  solutions. 

For  near  term  application,  a  viscous  ramp  procedure  is  suggested 

wherein  the  viscous  effects  are  phenomenologically  represented  by 
unsteady  shape  changes  that  can  then  be  handled  by  an  inviscid 

procedure.  On  a  longer  time  frame  a  more  sophisticated  lag  entrainment 
integral  method  is  suggested.  In  the  incorporation  of  the  viscous 

effects,  the  3D  flow  is  first  assumed  to  be  quasi-planar  in  a  suitable 
swept  framework;  and  in  the  latter,  the  viscous  flow  is  assumed  to  be 
quasi-steady  with  respect  to  the  instantaneous  inviscid  flow. 

The  coupling  of  the  lag  entrainment  method  with  the  transonic 
inviscid  procedure  will  be  a  challenge,  particularly  if  there  is  a 

stringent  limitation  on  the  added  computer  time. 


Finally,  one  must  not  forget  the  resonant  wind  tunnel  wall 
interference  that  can  have  a  major  effect  on  aeroelastic  phenomena.  Its 
absence  must  be  assured  before  test  data  can  be  validated.  The  solution 
to  the  problem  formulated  herein  will  play  a  vital  role  for  this  purpose. 
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APPENDIX 


The  problem  as  formulated  in  Section  III,  is  now  derived  in  more 
detail  and  in  several  representations. 

1.0  PHYSICAL,  DIMENSIONAL  RECTANGULAR  COORDINATES  AND  THE  FULL  POTENTIAL 

1.1  Coordinates  and  Variables 

X,Y,Z  cartesian  coordinates  with  X  in  the  freestream 
direction  and  Y  in  the  "spanwise"  direction  (ft) 

T  time  (sec) 

4>  velocity  potential  (ft^/sec) 

P  pressure  (lb/ft2) 

f  density  (slug/ft3) 

y  ratio  of  specific  heats 

a.  freestream  speed  of  sound 

U.,  freestream  velocity  (ft/sec) 

P.  freestream  pressure  ( 1 b/f t2 ) 

p  freestream  density  (slugs/ft3) 


1.2  Pressure  Relation 


JL  =  (±)r 

P.  \  9.  I 

b-  ^ 


1.3  Flow  Differential  Equation 


(?)T  +  (p^)z  +  <'?4v)y+  (tty*  =  0 


Subscripts  denote  partial  differentiation,  and  p  is  given  in 
Section  1.2. 


1.4  Shock  Jump  Conditions 


Equation  of  the  shock  surface:  S(X,Y,Z,T)  =  X  -  S°(Y,Z,T)  =  0 
Normal  vector  to  shock  surface:  Nt  =  (H  S^.+  I  S^+  J  Sy+  K  Sz)|Ns| 

l»/  ■  Sj  *  ♦  4  ♦  s| 


H,I,J,K  are  unit  vectors  in  T,X,Y,Z 
directions  respectively 


Shock  velocity  (Ua  ,VS  ,Wj. )  (ft/sec): 


sT+  sx  14+  srvs*  szk  =  o 


Jump  conditions: 

[M-Nt]=  0 


(In  terms  of  S°  =  X  -  S(Y,  Z,T),  in  the  above  equations  set  Sx=  1 
s;~V  S“  =  -5t) 


1.5  Wing  Boundary  Conditions 

Equation  of  wing  surface:  B(X,Y,Z,T)  =  £- £°(X,Y,T)  =  0 

Wing  Tangency  Condition: 

=  0 
or 

s;+  #XB‘+  ^2°  -  <b2  =  o 

Wing  Kutta  Condition: 

P*flTC.Y,B',Tj  =  P7x„  ,Y,B",TJ 

(+  ■  upper  surface;  -  «  lower  surface;  TE  *  trailing  edge) 


1.6  Wake  Condition 


Equation  of  Wake:  W(X,Y,Z,T)  =  0 

No  Load  Condition:  [P]  =  P-  P  =  0  ^  ^  ^um*) 

Flow  Tangency  Condition:  (V^-  )  *  Nw=  0 

where  +  denotes  the  upper  or  lower  surface  of  the  trailing  wake  and 
Nw  is  the  normal  vector  of  the  wake  surface  HVt«-  IWx+  JVy+  K V2 

The  primary  difficulty  encountered  here  is  that  the  side  edges  of 
the  wake  roll  up  in  an  ever-tightening  "spiral"  in  the  downstream 

direction  making  it  difficult  to  fulfill  the  wake  conditions  at 

downstream  points.  Moreover  the  potential  approximation  here  ceases  to 
be  valid  as  diffusive  effects  enter  the  picture  in  the  tightening 

"spiral".  Here  either  a  "solid  core"  must  be  introduced  or  more 
drastically  the  wake  condition  must  be  imposed  in  a  quasi-planar 

fashion.  Such  downstream  oversimplified  modeling,  fortunately  does  not 
significantly  distort  the  flow  over  the  wing. 

1.7  Outer  Conditions  (Approximate) 

Far  upstream  boundary:  =  LLX 

Far  lateral  boundaries: 

<£>  =  0  (Solid  walls) 

n 

Far  downstream  boundary  (Quasi-planar  wake) 

u,<Vn.)=  o 
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1.8  Wind  Tunnel  Wall  Conditions  (Highly  Idealized) 


On  side  walls: 


On  upper  or  lower  walls: 

i  flX 


1?4-*o ;  solid  wall 
-r^oo  ,  k=o:  free  jet 
TV*00  :  slotted  wall 
K0=D:  perforated  wall 


2.0  PHYSICAL,  NON-DIMENSIONAL  RECTANGULAR  COORDINATES  AND  THE 
DISTURBANCE  POTENTIAL 


2.1  Coordinates  and  Variables 

x  =  X/c  y  =  Y/c  z  =  ^/c 

t  =  (UL/c;  T 

=  iOc.Y.Z.Tj/dic;-  x 

(all  lower  case  variables  are  non-dimensional  and  c  is  a  reference  chord.) 


2.2  Pressure  Relation  (through  second  order  in  <f) 

C„  =  ,  =  -2(%+ft)- 

2 


: 

* 

* 

| 

i 


2.3  Differential  Equation  (To  Second  Order) 


where 

i2=  {'-H^VV} 
f}  =  {<-  «»V<p/n)}  ^  . 

2.4  Shock  Jump  Conditions  (To  Second  Order) 

J 

Equation  of  shock  surface:  $(*  >y  ,1  ,t)  =  X  -  S’Y.y.Z.'t)  -  0 

Shock  jump  conditions; 

[f.l  st  +  tf.]  s,  +  Mils,  +  [-fsJ 

I'f.)  s.  -  l*tK 

[%Ut  -[%]*, 

[9,i  st 

Shock  velocity  (Us, Vs, Vs): 

St  +  us  5,  -*  vs  sy  ■+  w5  52  =  0 
Here  Us  *  US/{J„  }  e+c.. 


=  0 
=  0 
-  0 

=  0 
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2.5  Wing  Boundary  Condition  (To  First  Order) 


Equation  of  wing:  2-  b(V,y,t)-  0 

Tangency  condition: 

b‘+  b°  -  =  0 

*t  X  * 

Kutta  condition: 

V'rE  ,y.0*,t)+  %(*rt  .V  ,o",  t)  = 

hfat .y.o'.t)-*  \itn  ,y,o~t). 

(0+  =  upper  surface;  0  =  lower  surface;  TE  =  trailing  edge). 

2.6  Wake  Condition  (Planar,  First  Order) 

No  load  condition: 

Tangency  condition: 

%(*■], o*,t)  =  %(*•*.<>'<  t) 
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2.7  Outer  Conditions 


Far  upstream  boundary:  f  -  0 

Far  lateral  boundaries: 


-  0 


Far  downstream  boundary: 

V  %  =  0 

2.8  Wind  Tunnel  Wall  Conditions  C  K  -  U„c  K  >  v  =  ) 

Sidewalls:  -2(%+<ft)  +  K  "*■  %  ~  ® 

Upper  or  lower  walls: 

-2(<PS%)+  k  r  0 

3.0  SCALED,  NUN-DIMENSIONAL  RECTANGULAR  COORDINATES  AND  THE  DISTURBANCE 
POTENTIAL. 

3.1  Scaling 

x  =  x  y  =  y/v  ?=  z/> 
t  =  tAc 

fyd.y.z,*)  =  6  ?(*,  J.z  .*  j 
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3.2  Pressure  Relation 


=  -2<=(vn)-£j{^v  "‘V  *r‘Vj* 


3.3  Differential  Equation 


(?.){  +  (f  )jr  (4V 


(fj?  - 


2  =  -  61  »!(%,+ 2vfy-  j-(‘l  M„z  2?$} 

+  {l-  (2-r) M*  }  ?»-*  +  -*  K  V 

£  =  fcO-Hi^y  -  -k^  [-^-r)X2Mw2  ^ 

+  2t:{i-  j^.  ^ 

+  (*+»)  V 

+  v2$-  +  K V  ]  *  •'* 

?3  =  +"'. 
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3.4  Shock  Jump  Conditions 


Equation  of  shock  surface: 


S(7.  Z  ,  t)  =  0 


Jump  conditions: 


ffJSf  +  V  LfJ  %  ♦  RJ  5,  =  o 


[  %]  S{  -  [%]•$* 


=  0 


=  0 


Shock  velocity 


L9s]Se 


-l%]$  =  0. 


~t  u  i  „  ■+ 

'"’t  ^  * 


~  ~  ~  7  r\ 

vi  Ss  +  Ws  st  =  0 


$.*  Sll,./  ^Us  ;  **■  ■l''5 :  5>=^Wi 


3.5  Wing  Boundary  Condition 


Equation  of  wing:  b°(xJ  i  )  -  0 


Tangency  condition: 


Kutta  condition: 


'v.  C  ~  . 

%  =  -2-  (  W  ♦  T  b*  ) 
i  6  jx  t  1  ' 


%lln,5,0',t)+  ziOn.S.o'.t) 


'  ' 

-  >»  .i/Tvy-.. 


3.6  Wake  Condition 


No  load  condition: 

\(x  j  t0*  +  ^  y,  0\t) 

=  %(*.?, *  fy*,y,o]t) 

Tangency  condition: 

3.7  Outer  Conditions 

/v- 

Far  upstream  boundary:  ~  0 

Far  lateral  boundary: 

ft,  *  0 

Far  downstream  boundary: 

%  +  t  %  -  0 

3.8  Wind  Tunnel  Wall  Conditions 

Side  walls:  ^  .  .  .  ~  / 

-2 ( -f  t fc.  ~  ^  y 

*  t  yy  R 

Upper  or  lower  walls: 

-  j(v  J  -£/***•  =  0 
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4.0  SPECIAL  CASE:  TRANSONIC  SCALING 


4.1  Scaling 


\)=  <rS* 


r-s*  -  z"/3 


i 
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6  -  si  ;  o-m*1)  -  k  <5 3  . 


4.2  Pressure  Relation 


Cp  =  -  % 


4.3  Differential  Equation 


(i)t~  +  (f')T  +  +  =  0 


l  - 


i,=  $T  J-  )c9>7-^M^r.0?/j 

V-(2-t,px‘  VJ 

f  -VsX*  ?,%•  -fie'S*  fj?t  *  -■ 


£  =  Sf  ?r-  -  fill  •  . 
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x'jafc  «v«*  . 


4.4  Shock  Jump  Conditions 


Identical  to  Section  3.4  with  ,  f,  ,  f,  ( ^ j 

obtained  from  Section  4.3. 

' 

Wing  Boundary  Condition 

Equation  of  wing:  Z  -  &*  b  ~  ® 

Tangency  condition:  ~  r4  To 

=  bj  -»  pS3  \ 

Kutta  condition: 

5  1 

+  P  S'3{tLK',y,o*,i)- 

A-  —  A<*  l) 

4.6  Wake  Condition 


4.7  Outer  Conditions 


Far  upstream  conditions:  *P  -  C) 

Far  lateral  condition: 


<f~  =  O 
n 


Far  downstream  boundary:  -i  ^ 

V  p*1  <f.  =  o 

4.8  Wind  Tunnel  Wall  Conditions 


Side  walls: 


0 


Upper  or  low  walls: 


-a?* «  V  0 


4.9  Comments 

If  in  the  above  equations,  we  take  the  limit  as  8  -»  0,  we  obtain 
the  low  unsteadiness  case  defined  in  Section  III  Subsection  7.0,  with  the 
spanwise  effects  represented  only  by  the  term.  There  is  no  scaling 
of  the  time  to  obtain  the  high  unsteadiness  case  with  the  <Pt  term  in 
this  limiting  process. 


The  corresponding  limiting  shock  jump  conditions  are  those  labeled 
Guderley-von  Karman  in  Figure  1;  and  as  seen  previously,  they  are 
completely  inadequate  for  swept  shocks  being  satisfactory  only  for  shocks 
essentially  normal  to  the  freestream. 

If  we  examine  the  effect  of  the  transonic  scaling  on  the  wing 
planform  shape  distortion  in  the  limit  8-*  0,  an  untapered,  unswept 
wing  planform  of  infinite  aspect  ratio  is  obtained,  so  that  the  limiting 
shock  jump  conditions  are  not  inconsistent  with  the  distorted  planform. 

It  is  amply  clear  that  no  scaling  exists  for  which  a  physically 
meaningful  limiting  process  would  yield  the  basic  flow  equation,  Eq.  (4) 
suitable  for  swept  wings. 
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APPENDIX  B 
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where  g  is  the  determinant  of  the  metric  tensor  of  the  sheared  transformation 
(B2)  and  is  given  by  g  =  c  (y)  (see  Lass,  p.  279);  and  f0  ,i,  )ii  ,f3  are  the 
transformed  components  of  the  contravariant  mass  flux  vector  given  by  (see 
Lass,  Eq.  491) 

l  -  i 


(B4) 


where  4y  from  Eq.  (B2)  is 


primes  denoting  differentiation  with  respect  to  y. 


Finally  in  Eq.  (B4)  we  must  transform  the  derivatives  arising  on  the 
right  side  with  the  result 

1=  *  V 

^  f VfcW-i?,1-  ir'H 

l  =  »tf  \ 


Thus  Eq.  (B5)  together  with  Eq.  (B3)  with  y/g7*  c  (y)  constitute  the 
transformed  equation  in  conservation  form.  It  is  the  proper  conservation  form 
with  the  components  given  in  Eqs.  (B4)  and  (B5)  representing  components  of  the 
mass  flux  vector. 
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M«o  *  0.84  —o—o-  Experiments 

01  *  6°  -  Calculations  (Ref.  2) 

ONERA  M-6  WING  (F™e  mesh) 

H  ~  0.45  - Calculations  (Ref.  2) 

(Coarser  Mesh) 


Figure  2.  Forward  Shock  Capture  with  NLR  Equation 


Figure  5.  Far-Field  Condition  -  Hyberbolic  Case 


$(x0,  y,  t  +  At)  *  ♦(xq  -  1/kAt,  yf  t) 


or 


k$t  +  ♦x  *  0 


Figure  6.  Far  Downstream  Boundary  Condition  (Trefftz  Plane) 


Figure  9.  Viscous  Interactions  -  Supercritical  Airfoil 


Figure  11.  Green's  Lag  Entrainment  Method  in  the  Case  of  Seddon's  Flow 
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